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ABSTRACT 

The observed orbital parameters of the recently discovered binary pulsar J0737- 
3039 are used to constrain the progenitor system. In particular, the small observed 
eccentricity and the small inferred peculiar velocity imply that during the formation 
of the second NS in the system, a very small mass AAf < 0.15Mq, was ejected. A 
progenitor more massive than IAMq is unlikely and a progenitor more massive than 
2.3Mq is practically ruled out. We therefore argue that the companion removed most 
of the progenitor envelope prior to its collapse. The second collapse may have been a 
SN of Type lb, or more probably, of Type Ic. Alternatively, it could have been a bare 
collapse, an option which is kinematically favored. Future constraints on the peculiar 
motion could solidify this constraint on the mass loss. A future determination of the 
pulsars' spin vectors and proper motion could allow us the complete determination of 
the conditions during the second collapse. 

1. Introduction 

The discovery of the new binary pulsar J0737-3039 (Burgay et al. 2003; Lyne et al. 2004) 
opened a new ideal general relativistic laboratory. The eclipses of J0737-3039A beyond J0737- 
3039B provide a superb way to explore pulsar magnetospheres. The relatively short life time of 
this system to gravitational radiation emission has lead to a revision of the binary merger rate in 
the galaxy (Burgay et al. 2003; Kalogera et al. 2004). We point out in this letter that the orbital 
parameters of this system, and in particular its low eccentricity and its low peculiar velocity pose 
strong limits on the origin of this binary system. We argue that the progenitor of the second 
neutron star must have been very light and that the formation of this neutron star must have 
involved only a minimal mass loss. 

The mass loss during the second explosion could easily disrupt the system. A mass loss prior 
to the second core collapse (usually due to a common envelope phase) or a kick velocity are usually 
invoked as ways to prevent the disruption of the binary during the second Supernova. Indeed there 
is evidence for kicks given to neutron stars at birth (see van den Heuvel & van Paradijs 1997; 
Burrows 1998). However, a kick velocity could either keep the ellipticity small or reduce the CM 
motion. That is to say, a kick which keeps the ellipticity small will increase the CM motion and 
vice versa. Therefore, the prior mass loss must have been very significant in this system. 
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The analysis described here is similar in spirit to that of Wex, Kalogera, &: Kramer (2000) who 
calculated the pre-SN characteristics of the binary pulsar system B1913+16 given the kinematics 
of the system. Nevertheless, several difference between the systems imply that both the analysis 
and the type of conclusions will be different. For example, we cannot place a limit on the direction 
of the natal kick relative to the pulsar spin direction. We can however place a very strong limit on 
the progenitor mass. 

On the other hand, the analysis here is different from two previous analysis of J0737-3039 
by Dewi k, van den Heuvel (2004) and Willems k. Kalogera (2004). In these analyses, the kick 
velocity of the neutron star was constrained given the center of mass kinematics of the system and 
evolutionary constraints. In particular, it was assumed that the pre-SN star was more massive than 
2.IM0 (Willems & Kalogera 2004) or 2.3M0 (Dewi & van den Heuvel 2004), since evolutionary 
scenarios do not give NS's forming from less massive pre-SN Helium stars in tight binaries. Here, 
we do not assume any evolutionary constraint, or for that matter, any assumption on the type of 
the progenitor. That is to say, we use here only kinematic constraints. Moreover, we also consider 
the kick velocity of the system as a whole, based on a probability analysis for finding the system 
as close to the galactic plane as it is. 

A minor difference between the analysis here and that of Dewi & van den Heuvel (2004) and 
Willems &; Kalogera (2004) arises from the fact it is now known that the companion of J0737-3039A 
is a pulsar as well, with a shorter characteristic spin down age. As a consequence, the last collapse 
in the system was more recent than assumed before. Nevertheless, the slow evolution of the orbital 
separation and eccentricity implies that this will only give rise to a small quantitative correction. 

We discuss the properties of the binary pulsar in §2. In §3 we estimate the probability that the 
binary system obtained different natal kicks. In §4 we consider the evolution of the system and the 
dynamics of formation of the second neutron star. We discuss the conclusions and possible caveats 
in §5. 

2. Pulsars' Properties 

The binary system J0737-3039 is composed of two pulsars denoted A and B. The system 
was discovered during a pulsar search carried out using a multibeam receiver at the Parkes 64-m 
radio telescope in New South Whales. At first only J0737-3039A was detected (Burgay ct al. 2003). 
Further observations revealed that the companion in this system is also a pulsar J0737-3039B (Lyne 
et al. 2004). We begin by summarizing the relevant parameters of this system. 

• Orbital Parameters: The orbital period is 2.454 hours. The separation of the pulsars today is 
8.8 X lO^'^cm and the eccentricity is 0.087779(5). This implies (see Fig. 2) that the system will 
merge within 83 Myr due to emission of gravitational radiation. The system is seen almost 
edge on and its inclination is 87°. 
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• Pulsar life times: The pulsars' periods, Pa,b are: 22.7 msec and 2.773 sec respectively. 

Pa,b and the time derivatives Pa,b provide estimates for the life times of the pulsars with 
~ 210 Myr and ~ 50 Myr. These observations also provide estimates for the magnetic 
fields of the two binaries: Ba ~ 6.3 x lO^G and Bb ^ 1.6 x IO^^q^ 

• Position: The system is ~ 600 pc away from Earth. Its Galactic coordinates are: L=245.2357 
and B= -4.5048, which implies that it is within 50 pc from the Galactic plane. 

• Masses: The pulsars's masses are almost equal with: mA = 1.377(5)M0 andMg = 1.25O(5)M0. 

We assume in the following that the present dual pulsar system was born about 50 Myrs ago 
when pulsar B was formed (as inferred from its lifetime). We infer that prior to that, pulsar A was 
spun up by accretion from its companion and that this spin up process stopped 200 Myr ago. 

3. Peculiar Motion 

One of the important observations which can be used to constrain the system its observed 
location — about 50 pc from the galactic plane. This is obtained given the dispersion measure 
distance of ~ 0.6 kpc and its galactic b = —4.5° latitude (Lyne et al. 2004). This implies that it 
is unlikely that the system attained a large peculiar motion during the core collapse. If it had, it 
would have left the galactic plane unless, by chance, it was kicked within the galactic plane. A 
large velocity would have made the system unbound from the Galactic potential. It would have 
then travelled, within its lifetime a distance of 5kpc(u/100km/sec)(t/50Myr). 

If the system was vertically bound then the motion would have been periodic. The vertical 
angular speed is given by 

nl = 4TrGpo. (1) 

A typical value for the density in the disk, po, is O.25M0/pc'^ (Bahcall et al. 1982). This gives an 
orbital period of Pz ~ 50 Myr. The typical velocity for an object observed at Zobs is 14 ~ ^T^Zobg/Pz- 
A distance of 50 pc would imply then that the vertical velocity is of the order of 6 km/sec. 

To quantify the probability, we perform a Monte Carlo simulation. We assume that the 
progenitor distribution has a Gaussian distribution in the amplitude Az^ of the vertical oscillation, 
with a width = 25, 50 or 100 pc. At the moment of core collapse, we assumed it has a phase cpo 
within its vertical motion, such that: 

2ttAzo , , 

= AzfiSiatpQ , Vz,o = — ^-^cos(/7o- (2) 

We have no particular priors on the kick velocity, we therefore assume the center of mass kick 
velocity has a log- normal distribution between 1 and 1000 km/s, and that the direction of the kick 
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is random. Given a random center of mass kick velocity Vcmi we then calculate the new vertical 
velocity, and vertical amplitude: 




Vz,l = Vz,0 + Vcm,z , Az,l = Y ^0 + ( ) 

The final vertical location of the binary will be given by: 

zi = Az^isinifi, (4) 

where we assume that tpi is random. If we could have been certain about the pulsar's age as well 
as the periodicity Pz, we could have calculated it from the kick and tpo, but this is not the case. 

We repeat this procedure 10^ times, and collect all the events which leave a binary located at 
zi = Zobs = 50 ± 5 pc. We then plot in Fig. 1 the distribution function P{vcm > Vcut)- That is, we 
plot the probability that the given system had a CM kick larger than Vcut- We find that, as expected, 
^cm ^ 15 km/s at 68% confidence, Vcm ^ 150 km/sec, at 95% confidence, and Vcm ^ 500 km/s at 
99% confidence. This velocity is small compared with the internal circular velocity of the system 
~ 600 km/s. This will prove to be important. 

The limit on the vertical velocity and on one of the perpendicular (to the line of sight) com- 
ponents of the velocity within the Galactic disk can be found within a few years, as estimates of 
the proper motion of the system will be obtained. Note that at a distance of 600 pc, a velocity of 
100 km/sec implies a proper motion of 0.036"/yr. 



4. Evolution 

At present the orbital evolution of the two pulsar's system is determined by emission of grav- 
itational radiation. This causes a spiraling in of the binary orbit. The orbital evolution is well 
known and can be integrated forwards and backwards in time according to the formulae given, for 
example, by Shapiro and Teukolsy (1983). Fig. 2 depicts the evolution of the orbit due to gravi- 
tational radiation emission. An interesting feature of the evolution is that the eccentricity decays 
with a faster rate than the orbital separation. Thus an eccentric orbit becomes, first circular and 
then it decays. 

Fig. 2 shows that the system will merge within 85 Myrs. This makes it the shortest lived known 
binary pulsar system. This short life time has led to a revision upwards of the estimated rate of 
binary neutron star mergers in the Galaxy by an order of magnitude (Burgay et al. 2003; Kalogera 
et al. 2004). Fig. 2 also shows the integration backward in time. This enables us to determine 
the parameters of the system when it was born. We find that 50 Myr ago, when the pulsar B was 
born, the ellipticity was shghtly larger (cj ~ 0.11) and the semi-major axis was ~ lO^-'^cm. Thus we 
conclude that the system was born with a low ellipticity and the present low ellipticity is not due 
to gravitational radiation decay of the ellipticity. Instead, it resulted from the orbital parameters 



-5- 



of the system after the second core collapse. Note that the values are not that different at 200 Myr. 
At that time the ellipticity was 0.14 and the separation was 1.2 x 10^^ cm. 

We come now to our main consideration. We ask how can the system form with both a low 
eccentricity and a small velocity? We show that this implies that the mass loss during the second 
core collapse was minimal. Wc assume that prior to the second core collapse, tidal interaction 
between the progenitor and the neutron star has led to a circular orbit. We consider the influence 
of mass ejection during this core collapse event on the orbital motion. We begin with a simple 
analytic toy model in which we ignore the eccentricity and demand that the binary neutron star 
system moves in circular orbit. 



4.1. A Toy Model 

Consider two masses rrii and m2 in a circular motion with a separation A^, before the NS 
formation. The star i explodes and sheds some of its mass. The remaining mass is m/. Within the 
context of J0737-3039 mf = niB and m2 = tua- We assume that the explosion is non-spherical and 
as a result the star attains a kick velocity Vk in opposite direction to its absolute velocity. Prior to 
the explosion, the relative velocity of rui relative to m2 is: 

2 G{m2 + rui) 
'"c,i = J. ■ (5) 

In order for the new orbit to be circular, the relative velocity of m/ relative to 7712 should be: 

^,^^G(™^ (6) 

This leads to a condition for Vk'- 



Y Ai y A-i 

The mass loss and the kick lead to a center of mass velocity: 



G(m2 + rrii) / G{m2 + ruf) 

-Vk = \ ■ (7) 



{nif + 1112) {rrii + m2) '^'^ mf + m2 

Defining: q = m2/m/; / = mf/rrii and ( = Vcm/vcj we get: 

C = I \- (9) 

from which it is self evident that with qK. 1, in order that C should be small (small CM velocity) 
/ must be close to unity. The results can be understood in the following way. A spherical mass 
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ejection gives a kick to the CM of the remaining system direction as the kick given to the CM. The 
kick is of the order of the circular velocity times the fraction of the ejected mass divided by the 
total remaining mass. Now, in order that the system remains in a circular orbit the ejected mass 
should be ejected non spherically in such a way that it gives a kick to m/, lowering its velocity 
relative to m2. This kick should of the order of the circular velocity of the initial system if the mass 
loss is significant. However, this adds a kick to the CM of the remaining system. This additional 
kick is in the same direction as the first one, leading essentially to a CM velocity which is of the 
order of the initial circular velocity. 



4.2. The realistic eccentric case 

We consider the same system again, with two masses mj and m2 moving on a circular orbit 
before the core collapse. Now the system attains an elliptic orbit after the SN. A will now be 
the sum of the orbits semi-major axes, or separation in case the pre-collapse circular orbit. The 
analysis uses the same principles as Hills (1983); Kalogera (1996); Fryer &; Kalogera (1997); Wex, 
Kalogera, & Kramer (2000). Here, we will obtain equations for the unknown initial progenitor mass 
rrii, the unknown initial separation Ai and the kick velocity Vk, given the current masses, orbital 
configuration and binary CM velocity. 

The relative velocity of rrii before and after the SN are y^rei,i and VrelJ = '^rel,i + Vk, where 
Vfc is the kick velocity. The condition of circular motion before the core collapse gives: f^^^ ■ = Vc,i, 
where Vc,i is given by Eq. 5. Defining 6 as the angle between and Vr^i^i and the orthogonal angle 4> 
to be zero if Vk lies within the plane defined by the momentary separation between the two masses, 
r, and Vrei,i, and pointing outwards from r, we have, 

VrelJ = \J {Vrel,i + Vk COS 9^ + {vk sin 9^ . (10) 

The energy of the system after the core collapse is: 

E = -g"^^^^ + [{vrei,i + Vk cos Of + {vk sin ef] . (11) 

where /x = mfm2/{mf + m2) is the reduced mass (after the core collapse) and A^ is the separation 
just prior to the core collapse (because before the collapse, the orbit was circular). The angular 
momentum and the eccentricity e after the core collapse are: 

Lf = AiH^ {vrel,i + Vk cos 9)'^ + {vk sin 9 sin cf))"^ (12) 

and 

^ "^^G2(m/ + m2)/x3 ^^^^ 



The new separation Aj^ satisfies: 



A, = -^<!^I±^ (14) 
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Finally, the velocity of the center of mass, Vcm after the core collapse is given by: 

rrii — mf\ f m2 \ mj 



mf + m2j \mi + m2 
Using our definitions of angles, we have: 



nij — in 



f 



ruf + m2 



1712 



rui + m2 



ruf 



ruf 



^rel,i + 



Vk COS 9 



ruf + m2 



+ 



ruf 



Vk sm ( 



(15) 



(16) 



We define an equivalent circular velocity for the binary after the core collapse: 

2 _ G{m.f + 777,2) 



A 



f 



(17) 



This is a fiducial quantity that we will use as a measure of the various velocities in the system. 
Defining k = Vk/vcj, and a = Af/Ai, and using the previous definitions of /, q, and we obtain 
three equations for the eccentricity, the semi-major axis and the CM velocity: 



«V'(l + g) 



^^—2 (1 + g) _ a (/ (2 + g) - 1) + 2V^fk,/T+^^U+~q cos 9^ 



x(^a{l + fq) + 2V^fk^yl + q^J+q cos 9 + fk^ (1 + q) cos 6^+ 
fk'^{l + q) sm9^sm{4>f^ , 



a (/ (2 + g) - 1) - ifk^ (1 + q)) - 2^fkVT+qJ^ + gcos 



and 



«(/ - 1) + (1 + 5) (1 + fq) + 2^ if- 1) fkq,/TT^^j +qcos9 

f{l + qfil + fq) 



(18) 



(19) 



(20) 



For a given e, 6 and (/) a solution to the three equations can be found numerically for f,k, 
and a. Then, the 6 and (p phase space is searched to find the minimum value of / possible. This 
gives the maximal mass that can be ejected for a given CM velocity. 

Fig. 3 depicts, for various initial eccentricities, the minimal / value as a function of The 
integration backwards of the orbit to the time of the core collapse, about 50 Myr ago shows that 
just after the core collapse the orbital separation was ~ 10^^ cm and the eccentricity, was around 
0.11. The corresponding circular velocity was 590 km/sec. Given that the current CM motion is 
most probably of order 15 km/sec, C < 0.025 is a reasonable upper limit. Given the eccentricity 
0.1, this value of C implies that / > 0.88 or that the mass loss during the formation of the second 
neutron star formation was probably of order ~ O.ISMq. For < 0.25 (our 95% confidence limit), 
we find that fmin = 0.55, thus, the star shed less than 1.0 Mq at 95% confidence during the core 
collapse event. The progenitor was most probably less massive than 2.3Mq, with rrii ~ IAMq 
being more likely. 
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5. Conclusions and Implications 

The binary B1913+16 was studied for almost 3 decades since its discovery (Hulse & Taylor 
1975). The ample data could then be used by Wex, Kalogera, & Kramer (2000) to place very 
interesting limits on this system. For example, it was found that the natal kick must have been 
directed almost perpendicular to the spin axis of the neutron star progenitor. On the other hand, 
the rather large post-SN eccentricity implies that the pre-collapse orbital separation (1.8 — 4.6 Rq) 
or the mass of the progenitor (4 — 32 Mq) cannot be tightly constrained. 

The system J0737-3039 is different from previously detected binary pulsar systems in that it 
has a very low eccentricity, and it was small also after the collapse. Thus, the orbital separation, 
for example, can be constrained to have been Ai = 1.0 zt 0.1 x lO^^cm = 1.45 it 0.15 Rq during 
the collapse. As the system J0737-3039 was only recently discovered, the proper motion could 
not be detected yet. Therefore, in order to get meaningful constraints on the kinematics of the 
SN explosion, Dewi &; van den Heuvel (2004) and Willems & Kalogera (2004) constrained in their 
analyses the progenitor mass to be consistent with the evolution of He stars in tight binaries. 

Statistically, however, it is possible to place limits on the center of mass motion of the system, 
at different confidence limits, using the fact that it is within about 50 pc from the galactic plane. In 
particular, Vcm ~ 15 km/s and it is less than 150 km/s at about 95% confidence. This implies that 
the ejected mass during the collapse is likely to be less than ~ 0.15 Mq, while at 95% confidence, 
it is less than ~ 1.0 Mq. 

Since current evolutionary models for Helium stars require a minimum mass of 2.1 to 2.8M0 to 
form a NS (Tauris & van den Heuvel 2004, and references therein), there seems to be two reasonable 
possibilities. 

In the first possibility, J0737-3039B originated in a SN by a progenitor Helium star with 
Mq > 2.1Mq. This case requires the systemic kick to have been large but within the galactic 
plane. This could have occurred with a low probability of < 5%. Since He stars which lose 
mass only through winds generally have somewhat higher masses ranging between 2.3 to 3.5 Mq 
(Woosley, Langcr, & Weaver 1995), a possible interpretation is that the progenitor star lost most 
of its envelope through interaction with its companion. This is also reasonable considering the 
small separation. Since most of the lost mass could not have accreted onto the companion, it was 
probably lost through a common envelope phase, at which point the companion J0737-3039A was 
spun up and its magnetic field was suppressed by accretion. This accretion phase ended around 
200 Myr ago. 

Because so little of the star was left after the mass loss, it is possible that even the Helium 
mantle was removed, such that the resulting supernova, if it were a supernova, was of Type Ic. 

A second interpretation of the above results follows the more probable mass loss. We expect 
the mass loss in the cataclysmic event to have been < O.ISMq, given the small systemic velocity. 
Moreover, we note that just the u burst from a core collapse should reduce the mass by the same 
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order, with a typical value being ^ IQ^^erg/c? ~ O.OSMq. In other words, the collapse may have 

been of a naked core. Such a collapse does not eject any appreciable mass, but it sheds enough 
"mass-energy" in the form of u's to explain the small eccentricity. In such a case, no kick velocity is 
expected, and we can determine the v burst energy from the inferred eccentricity after the collapse. 
Namely, E^, « e(l + q)/[{l + e(l + q)]mBC^ = Q.22S>t Q^Q^r,MQC^ k, 4.2 x lO^^g^^ (assuming the 
collapse took place 50^25 ^y^" ago). If some mass is ejected, E^, will of course be smaller but the 
total mass-energy lost will be the same. Thus, under the kinematically more probably option, this 
collapse did not proceed along the common core collapse scenario. Instead the WD progenitor, with 
a mass of Mg + Eyjc? = 1.478j]QQ^5MQ was supported by thermal pressure as well as degeneracy 
and that the collapse was induced due to the loss of this thermal pressure as the progenitor was 
cooling. It could have also ensued other internal transitions such as the process of slow e-capture 
which reduces the Chandrasekhar limit (Finzi k. Wolf 1967). 

At present we cannot, however, rule out the kinematically less probable possibility that the 
system does have a large CM velocity and that this velocity is within the Galactic plane. This pos- 
sibility implies a He star progenitor for J0737-3039B . The kinematically more favorable possibility 
implies that the pulsar is the remnant of the naked collapse of a WD. 

The future detection within a few years of a proper motion could constrain the center of mass 
velocity much more, after an additional velocity component would be determined. A constrained 
peculiar motion will probably solidify this extreme case of a low mass collapse and resolve the 
puzzle on the origin of the object. A useful measurement in this respect, is of the inclination of 
the pulsar orbit relative to the galactic plane. This could be achieved through measurements of the 
polarization behavior and of the relativistic spin precession. 

The authors wish to thank Ramesh Narayan and members of the Israeli center for High Energy 
Astrophysics, D. Eichler, M. Milgrom, E. Waxman and V. Usov for helpful comments. The research 
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Fig. 1. — The probability that the binary attained a peculiar motion given by Vcm after the core collapse, that is 
larger than a cutoff speed Vcut, given that today, it is located roughly ~ 50pc from the galactic plane. The dashed, 
solid and dotted lines assume that az = 25,50, 100 km/s. The two latter cases are virtually indistinguishable. 
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Fig. 2. — The evolution of the binary "separation" A defined as the sunn of the semi-major axes of both pulsars 
and the orbital eccentricity e due to the emission of gravitational waves as a function of time. This was achieved 
by integrating the orbital parameter evolution equations (Shapiro & Teukolsky 1983) both forward and backward 
from the current configuration with A = 8.8 x lO^^cm and e = 0.088. 
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Fig. 3. — The minimal fraction of mass / = mf/rm left after the core collapse in units of the initial mass rrii of 
the primary, plotted as a function of the kick velocity |vcm|, or in a dimensionless velocity ( = |vcm|/wcirc,/- The 
different graphs are calculated for different eccentricities, with e sa 0.1 being the preferred value at t ~ —50 Myr, 
the rough age based on the spin down age. ( is expected to be ^ 0.02. Otherwise, the binary would have been 
far from the galactic plane, unless by chance it was kicked very close to the plane. 



